arXiv: gr-qc/0512051 v3 15 Feb 2006 


February 5, 2008 15:51 WSPC/INSTRUCTION FILE 


br4 


Bel-Robinson tensor and dominant energy property in 
the Bianchi type I Universe 

BIJAN SAHA*, VICTOR RIKHVITSKY 
Laboratory of Information Technologies 
Joint Institute for Nuclear Research 
141980 Dubna, Moscow region, Russia 

MIHAI VISINESCU t 
Department of Theoretical Physics 
National Institute for Physics and Nuclear Engineering 
Magurele, P.O.Box MG-6, RO-077125 Bucharest, Romania 


Within the framework of Bianchi type-I space-time we study the Bel-Robinson tensor 
and its impact on the evolution of the Universe. We use different definitions of the 
Bel-Robinson tensor existing in the literature and compare the results. Finally we inves¬ 
tigate the so called ’’dominant super-energy property” for the Bel-Robinson tensor as a 
generalization of the usual dominant energy condition for energy momentum tensors. 
Keywords: Bianchi type I model, super-energy tensors 
Pacs: 03.65.Pm and 04.20.Ha 


1. Introduction 

The lack of a well-posed definition of local energy-momentum tensor is the conse¬ 
quence of the Principle of Equivalence^, which lies at the heart of Einstein’s theory 
of general relativity. Nevertheless, quest for the local tensors describing the strength 
of gravitational held has long been going on. One of the first successful attempt to 
address this problem was taken by Bel 2 3 4 an d independently Robinson 13. In 
their works, in analogy with the electromagnetic energy-momentum tensor, they 
constructed a four-index tensor for the gravitational held in vacuum. The prop¬ 
erties of the now famous Bel-Robinson (BR) ten sor are similar to the traditional 
energy-momentum tensor and following Senovilla ® "^can be formulated as follows: 
(i) it possesses a positive-dehnite time-like component and a ’’causal” momentum 
vector; (ii) its divergence vanishes (in vacuum); (iii) the tensor is zero if and only 
if the curvature of the space-time vanishes; (iv) it has positivity property similar 
to the electromagnetic one; and some others. Construction of BR and the study 
of its p roperties were widely considered by a number of auth ors, e.g., Desher et. 
al. 12“ Teyssandier D21 Senovilla F Bergqvist Hlj Andersson Wingbrant CCS 
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Choquet-Bruhat it et. al.^^etc. 


It should be noted the authors of the papers mentioned above considered the 
BR and established its properties in general. On the other hand, in general relativ¬ 
ity there exists a number of interesting and widely studied models of space-time. 
Therefore, in our view it is interesting to consider the BR within the scope of some 
concrete metric. In a recent paper we studied the BR within the framework of 
Bianchi type I (BI) universe using two different definitions. The purpose of this 
paper is to extend that study for some other definitions and analyze the dominant 
energy property (DEP) and dominant super-energy property (DSEP) within this 
model. 

2. Bianchi I Universe: a brief description 

A Bianchi type-I (BI) universe is the straightforward generalization of the flat 
Robertson-Walker (RW) universe and is one of the simplest models of an anisotropic 
universe that describes a homogeneous and spatially flat universe. It has the agree¬ 
able property that near the singularity it behaves like a Kasner universe, even in 
the presence of matter, and consequently falls within the general analysis of the 
singularity given by Belinskii et al. Also in a universe filled with matter for 
p = £ e, <C < 1, it has been shown that any initial anisotropy in a BI universe 
quickly dies away and a BI universe eventually evolves into a Friedmann-RW (FRW) 
universe ^. Since the present-day universe is surprisingly isotropic, this feature of 
the BI universe makes it a prime candidate for studying the possible effects of an 
anisotropy in the early universe on present-day observations. In light of the im¬ 
portance mentioned above, several authors have studied BI universe from different 
aspects. 

A diagonal BI space-time is a spatially homogeneous space-time, which admits 
an Abelian group G 3 , acting on spacelike hypersurfaces, generated by the spacelike 
Killing vectors Xi = di, x 2 = <9 2 , and X3 = 83. In synchronous coordinates, the 


metric is ^ 


3 



( 2 . 1 ) 


If the three scale factors are equal (i.e., ai = 02 = 03), Eq. ( 12 . 11 ) describes 
an isotropic and spatially flat FRW universe. The BI universe has a different scale 
factor in each direction, thereby introducing an anisotropy to the system. Thus, a BI 
universe, being the straightforward generalization of the flat FRW universe, is one 
of the simplest models of an anisotropic universe that describes a homogeneous and 
spatially flat universe. When two of the metric functions are equal (e.g., 02 = < 23 ) 
the BI space-time is reduced to the important class of plane symme tric space-time 
(a special class of the locally rotational symmetric space-times - 0 ) , which admits 
a G 4 group of isometries acting multiply transitively on the spacelike hypersurfaces 
of homogeneity generated by the Killing vectors xi, x 2 , x 3 , and X 4 = x 2 d 3 — x 3 d 2 . 
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The BI has the agreeable property that near the singularity it behaves like a Kasner 
universe, given by 

ai(t) = a°t Pl , a-i (t) = afy ? 2 , 03 (f) = a 3 t P3 , 

with pj being the parameters of the BI space-time which measure the relative 
anisotropy between any two asymmetry axes and satisfy the constraints 

Pi+P2+P3 = l, pl + pl + p 2 3 = 1 . (2.2) 

As one sees, pi , P 2 and p 3 cannot be equal. Only two of them can be equal, and 
only in two special cases, namely, (0,0,1) and (—1/3, 2/3, 2/3). In all other cases 
Pi , P 2 and p 3 are different, moreover, one of them is negative, while the two others 
are positive. If it is supposed that p 3 < P 2 < P 3 , then their values are confined in 
the following intervals: 


—1/3 < _pi < 0, 0<p 2 <2/3, 2/3<p 3 <l. 


The solutions of the algebraic equations (l 2 ~ 2 l) can be presented as 

—p p{p +1) P+ 1 

P 1 2 i i T ' P^ 2 i i T 5 2 i i T* 

p+p +1 p +p + 1 p + p +1 

Thus instead of three, we have now one parameter p , which lies in the interval 

0 <p < 1 . 

Another particular parametrization can be given using an angle on the unit 
circle, since Eqs. (l 2 ~ 2 l) describe the intersection of a sphere with a plane in the 
parameter space (pi 1 p 2 1 p 3 ): 


1 + cost! + v^sim? 1 + cost? — \/3sini9 1 — 2cos$ 

Pl = -3-’ P2 = -3-’ P3 = -3-' 

Although 1 } ranges over the unit circle, the labeling of each p 3 is quite arbitrary. 
Thus the unit circle can be divided into six equal parts, each of which span 60°, and 
the choice of pj is unique within each section separately. For = 0, p\ = P 2 = | 
and p 3 = —3 while for 1 } = 7 r/ 3 , p\ = 1 and P 2 = P 3 = 0. 

For later convenience we list the Christoffcl symbol, scalar curvature, Ricci, 
Riemann and Weyl tensors for the BI space-time. The non-trivial Christoffel symbols 
are the following: 


r z _ 

1 zo — > 

CLi 


t 2, 3 , 


and the non-trivial components of Riemann tensors are 
R 0i _ 

i,j = 1,2,3, 


j-*ij _ Ctj 

01 = ij ~ 


(2.3) 


I 


l 
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Finally the nontrivial components of the Ricci tensors for the BI metric are 

* 0 ° = -E 

i= 1 



R\ = - 

di 

+ * 


+ -)1 

di 


.di 

di 

\dj 

afc/J 


iijjk — 1) 2,3, i ^ j k. 


and the scalar curvature is 


R=- 2 


O 

(E 


(M ai a 2 a 2 a 3 a 3 a\ 

CLi CL\ a 2 a 2 CI3 03 Oi 


)■ 


It is convenient to separate the Riemann tensor into a trace-free part and a 
’Ricci” part. This gives the Weyl tensor 


Cijkl — Rr-,kl 


1 


■ij kl 


(n-2) 


^QikRjl T Qjl Rik QjkRil QilRjk^ 


1 


(n — 1 )(n — 2) 


QikQjl QilQjk 


)r. 


This tensor has manifestly all the symmetries of the Riemann tensor; but con¬ 
trary to the Riemann tensor while it gives rise to Ricci tensor, the Weyl tensor gives 


9 ik C ijkl = 0. 


(2.4) 


A further distinction is that while the Riemann tensor can be defined in a man¬ 
ifold endowed only with a connection, the Weyl tensor can be defined only when 
a metric is also defined. In 4 dimensions the Riemann tensor has 20 distinct com¬ 
ponents, while the Weyl and the Ricci have 10 components each. The non-trivial 
components of the Weyl tensor for the BI space-time are 


( 2 c , 0.'y dik djdkdi dkdidj dkdidj I 2 djdj^dq^ 


Com = 7 - 

0 djdk 

Cjkjk = 7 “ { ‘2'(lj / djdk djdkdi dkdidj didjdk dkdidj T 2 djdkdi^, 

Odi 

i,j,k = 1,2,3, (2.5) 


From one easily finds the following relation: 

a 2 

Com = - T~2 Cjkjk, i,j, k = 1,2, 3, i / j ^ k. (2.6) 

a j a k 

Now having all the non-trivial components of Ricci and Riemann tensors, one 
can easily write the invariants of gravitational field which we need to study the 
space-time singularity. Moreover now we can construct the BR tensor that is defined 
differently by different authors. 
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3. Einstein equations and their solutions 

In this section we study the Einstein equation. In doing so let us first write the 
Einstein equation for the BI metric governing the evolution of the Universe. In 
presence of a cosmological constant A the Einstein equation has the form 


— + — + = kT} + A, 

(3.1a) 

a 2 a 3 a 2 a 3 

— + — + = kTo + A, 

(3.1b) 

«3 0 ,i 03 Oi 

di + d a + did a= 3 + 

(3.1c) 

a\ ct 2 a i a 2 

Or d 2 a 2 d 3 d 3 Oi 3 

- 1 - 1 - = Kir. + A. 

(3-ld) 

cl 1 a- 2 a 2 a 3 a 3 Oi 

Here over-dot means differentiation with respect to t and T is 
momentum tensor of the matter field which we choose in the form: 

the energy- 

T I1 = (e + p)u li u v -p6”, 

(3.2) 


where is the flow vector satisfying 

g )LV u^u v = 1 . 

Here e is the total energy density of a perfect fluid and/or dark energy density, 
while p is the corresponding pressure, p and e are related by an equation of state 
which will be studied below in detail. In a co-moving system of coordinates from 
El one finds 

T 0 ° = c, T{ = Tl = Tl = -p. (3.3) 

In view of El from 13.1 all - 13.1 dll one immediately obtains 

ai(t) = Di[T(t)\ 1/3 exp[A/ J \T(t')]~ l dt'], i= 1,2,3. 

Here Di and A) are some arbitrary constants obeying 

D^D 2 D 3 = 1 , Xi + X 2 + X 3 = 0, 
and r is a function of t defined to be 

t = a\a 2 a 3 . 

From 13.1 all - EH for r one find 

^ = y(--p)+ 3A - ( 3 -4) 

On the other hand the conservation law for the energy-momentum tensor gives 
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After a little manipulations from and 18.511 we find 

i 2 = 3 (ke + A )t 2 + Ci, 


(3.6) 


with Ci being an arbitrary constant. Let us now, in analogy with Hubble constant, 
define 


T 

r 


a b c 
- + T + - = ZH. 
a b c 


On account of cnt from (THU) one derives 


ke = 3H 2 — A — Ci/(3t 2 ). 


(3.7) 


(3.8) 


It should be noted that the energy density of the Universe is a positive quantity. 
It is believed that at the early stage of evolution when the volume scale r was close 
to zero, the energy density of the Universe was infinitely large. On the other hand 
with the expansion of the Universe, i.e., with the increase of t, the energy density 
e decreases and an infinitely large r corresponds to a £ close to zero. Say at some 
stage of evolution £ is too small to be ignored. In that case from CD9 follows 

3IL 2 — A —>■ 0. (3.9) 

As it is seen from (EH31) in this case A is essentially non-negative. We can also 
conclude from cm that in absence of a A term beginning from some value of 
r the evolution of the Universe comes stand-still, i.e., r becomes constant, since 
H becomes trivial, whereas in case of a positive A the process of evolution of the 
Universe never comes to a halt. Moreover it is believed that the presence of the dark 
energy (which can be explained with a positive A as well) results in the accelerated 
expansion of the Universe. As far as negative A is concerned, its presence imposes 
some restriction on £, namely, e can never be small enough to be ignored. In case 
of the perfect fluid given by p = (e there exists some upper limit for r as well (note 
that r is essentially nonnegat ive, i.e . bound from below). In our previous papers we 
came to the same conclusion — (with a positive A which in the present paper 
appears to be negative). A suitable choice of parameters in this case may give rise 
to an oscillatory mode of expansion, whereas in case of a Van der Waals fluid the 
highly nonlinear equation of state may result in an exponential expansion as well. 

Inserting and m into cu> one now finds 

H = — 2 ( 3 H 2 ~ A + ^2 + K P) = ~2 ( £ +P ^ ~ 3^2 ’ ( 3 ’ 10 ) 

In view of m, from Jinn . it follows that if the perfect fluid is given by a stiff 
matter where p = £, the corresponding solution does not depend on the constant 
Ci. 

Let us now go back to the Eq. (TOIl . It is in fact the first integral of and 
can be written as 


f = ± \JCi + 3 (ke + A )t 2 


l 


l 


(3.11) 
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On the other hand, rewriting 13.51) in the form 

£ f 
£ +p r ’ 

and taking into account that p is a function of £, one concludes that the right hand 
side of the Eq. £3> is a function of r only, i.e., 

Q /<* 

f = — [e — p)r + 3Ar = 5"(r). (3-12) 

From a mechanical point of view Eq. unit can be interpreted as an equation of 
motion of a single particle with unit mass under the force 5"(r). Then the following 
first integral exists 

f = V2[£-U(r)]. (3.13) 


Here £ can be viewed as energy and 'U(r) is the potential of the force 95 Comparing 
the Eqs. 13.111) and (I3T31) one finds £ = Ci/2 and 

U(t) = — ^(«£ + A)r 2 . (3-14) 

Let us finally write the solution to the Eq. <E3 in quadrature: 


dr 


\J Ci 3(k£ + A)t 2 


— t + to; 


(3.15) 


where the integration constant to can be taken to be zero, since it only gives a shift 
in time. The Eqs. E3> and 13.51) for p erfect fluid obeying different equations of 
state has been thoroughly studied by us ^ 


4. Bel-Robinson tensors 

BR tensor first appeared in the endless search for a covariant version of gravita¬ 
tional energy. In general relativity, the energetic content of an electromagnetic field 
propagating in a region free of charge is described by the well-known symmetric 
trace-less tensor 

Tf = F aX Ff - (4.1) 

where F a P is the electromagnetic field tensor. This tensor satisfies: 

J el;a — u 

as a consequence of Maxwell equations with j 1 ' = 0. The tensor T'*f enables us to 
define a local density of electromagnetic energy as measured by an observer moving 
with the unit 4-velocity u: 

w e \(u) = T^u a U 0 . 

It follows from ED that the energy density is positive definite for any time-like 
vector u. 
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Within the scope of general relativity, however, it is well known that the concept 
of local energy density is meaningless for a gravitational field. To overcome this 
difficulty led to introduce the notion of super-energy tensor constructed with the 
curvature tensor R^ap- The first example of such a tensor was exhibited by Bel 
El that was further generalized to the case of an arbitrary gravitational field 
Note that a similar tensor was also introduced by Robinson'^. This tensor is now 
commonly know as the BR tensor as well. Since we are going to compare some 
distinct definition of BR in this paper, before defining them let us see what kind of 
properties they should have. In general, the BR tensor has the following symmetry 
properties: 



(4.2a) 


(4.2b) 


(4.2c) 


The symmetry property leads to the fact that that in ro-dimensional case there are 
n(n+ l)[n(n + 1) + 2]/8 independent components of the BR tensor. In case of n = 4 
out of 256 components only 55 are linearly independent. 

In literature there are a few definitions of BR. Here we mention only three. 

I. By analogy with the tensor o which may be written as 

± l±v 1 fia. 1 i/ i n<y. ^ r i/ i 

the BR tensor is defined as^: 

Bfiuap = R P iT a Rpva0 + a * Rpvap- (4.3) 

Here the dual curvature is *R pu Xa . = (1/2) a pR a ^ Xa - It should be noted that 
this definition is adequate only in 4 dimensions and in vacuum. Otherwise this 
tensor cannot satisfy the DEP and therefore this expression should not be used 
in other dimensions or in non-Ricci-flat spacetimes. 

Using the definition of dual curvature, from (TOli we find 

= R P ^ a R P va@ + R? ° pRpvaa - -^9pvR a ^Rppar- (4.4) 

The properties 14. 2 all and 14.2hl) follow immediately from <EU thanks to the 
symmetry property of Riemann tensor. The property I4.2cl) is straightforward from 
(TOli . but for 14. 411 it requires 

gpuR a paT Rp P ar = g a pR/ aT R V par- (4.5) 

In view of 12.311 for the BR tensor in this case we obtain the following non-trivial 
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components: 


3 .. 2 

£>0000 — / o' 


B,, 


U a V 

= of + ctU % + % 

a j a l 


Bn 


CLj 


ak 


OiOi — CLi^i 4 } i r. 


Clo 


ak 


Bd q i q - 


= CliCLj + dididjdj—j , 
a k 

to -2 ■ 2 I a j 

J^ooti — a^ d ^ ^ 


Bujj — 


2*2 2 --2 2*2 0 >k 

= a i a j ~ a i a j ~ a i a j—i 
a k 


i,j,k= 1,2,3 

Inserting (EH1) into (TOli we obtain following additional relations: 


± 


2 r 


± 


2-i 


i,j,k = 1,2,3 i^j^k (4.6) 


Among the six constrains in (EH1 only three are linearly independent. After a little 
manipulations with them finally obtains the following relations between the metric 
functions: 

- = i,3,k= 1,2,3 (4.7) 

(Xj dk 

As one sees, in account of o the Einstein equation 14. la! - 15.1 dll leads to = 
T-l = T 2 2 = T 3 3 , which can be realized only when the source field satisfies the 
following equation of state: 

p = —£. (4.8) 

It is well known that only vacuum satisfies the state of equation given by lEDfll . 
Thus we see that if we are to define BR tensor given by o or foil we should 
deal with the Einstein equations with the source field given by a vacuum. 

II. The restriction that arises above is due to the fact that in defining the BR 
tensor we used the dual term with the duality operator act ing on the left pair only. 
To avoid this restrictions the BR tensor can be defined byl^" -^1 


2R M ;ya/3 = R p ^ aRpmyP + *R P £ a * Rpvaf) 

+ & * P g a R *pvaP + * R * P M Cr a *R*puop, (4.9) 

where the duality operator acts on the left or on the right pair of indices according 
to its position. Nowadays this is known as the Bel tensor and was introduced by 
Bel 2 in a slightly different form. 










February 5, 2008 15:51 WSPC/INSTRUCTION FILE br4 


10 B. Saha, V. Rikhvitsky, M. Visinescu 


From 14.911 one easily finds 


B 


yisa/3 


= R p „ 


yRpua(3 

par v 
1 1 1 


R P H pRpvcra 2 9p.v R o 

1 


r Rp, 


par 


c)9ot(3Rp P RvpcTT + QppvQoLfiR? ^Rparrj • 


Under the new definition the symmetry properties (TPal) . 14.2bl) and 14.2cll follow 
immediately, without any restriction to the metric functions. 

Let us now write the non-trivial components of the BR tensor for the BI metric. 
In view of (IQ) we now find 


Roooo — 


3 .. 2 

\ ^ ai_ 

aj 

1=1 1 


• 2 *2 
a{ 


• 2 *2 
0-2 ai 


a| a\ 


Baa — &i Boooo i 


BoiOi — CLiCLi 


Bijij CLiCLj 


Boon — 2 

_ 1 

Biijj — 2 


_J_ 

ij aj 


O'k Cik 
CLk &k 


*1'^") 9 

at 


9 ' „2 j 

| + a i 

(4 

, d l 

' 9 

a j o 

•2 



2 _2 -2 

a k t; 

,2„2 

i ^2 • 


a? al 


a2 J a l 


•2-2 -2 2 2-2 ^k (• 2 2 , 2 • 2\ , 2 2 u 'k 

a^ aj — dj aj — a^ dj — — ^ ®j ~F a^ dj j -\- a^ aj —^ 


i,j,k = 1,2,3 


But the BR tensor defined in this way is not trace-free and is not completely sym¬ 
metric. It is achieved if and only if the manifold is Ricci flat, i.e., Rij = 0. Since 
for the BI universe we have non-trivial components of Ricci tensor, we give an 
alternative definition of BR where it is totally symmetric and trace-free. 

III. Here we give another definition that gives rise to BR tensor, that is trace¬ 
less an d tota lly symmetric. It can be achieved by constructing BR by means of Weyl 
tensor 


B^a/3 = C P /J a a Cp Va p + *C P ° a * Cpva/3. (4-10) 

It can be shown that this BR is totally symmetric, i.e., 

Rijkl B('ijk:i )) 

Moreover, as one can easily find from 12.-41 1. the BR defined through Weyl tensor is 
trace-free, i.e., 

g :l B ljkl = 0. (4.11) 

Let us study this case in detail. Using the properties of Levi-Civita tensor we 
first rewrite iTHUl) in the form 

iW = C V + C p ; p C pV(Ta - l -g^C a O^CpprT. 
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After a little manipulation one easily finds the following nontrivial components of 
the BR tensor. 

-Soooo = (<? U Cioio) + (g 22 C 202 o) + (g 33 C 3030 ) } 
which on account of ESI gives 


60000 — 


(if 2323 ) + (i C3131 ) + (i^ 1212 ) ■ (4 - 12) 


1 


1 

a(a 2 


On the other hand for we find 

Bim = ^Oioioj + ^.9 22 Ci2i2^ + (<7 33 Ci3i 3 ) 
which in view of 12.011 can be rewritten as 

Bini = (^pi c2323 ) + (4 C3131 ) + (^ Cm2 ) 

( 9 2 ^ 2 323) + ( —o 2 ^313l) + ( 2 2 ^ 1212 ) 
\aiai / \aia1 / \a1ai / 






In view of and mnn symbolically we can write 


B 0000 — 


Ly k Cjkjk ) + Lla? Ckikl ) + (a?a? Cijij ) 


a i a j 


(4.13) 
= aiRoooo- 

(4.14a) 


= E 


3 .. 2 

af 


u, 2 




M 2 u 3 


h 2 3 


(a i a,2 

02 «3 

a 3 ai\ 

di a 2 , 

(a i 

, ®2 

0 d 3 

, a 2 \ 

1- 

1- 

H -) 



H - 

— 2 — 

+ ~2 

\d\ d2 

«2 a 3 

a 3 ai) 

ai a 2 

^ di 

a 2 

a 3 

a V 


a,2 a 3 , 

(02 


,o\ aj\ 

o 3 di , 

( d.3 

, Ol 

,02 

, at\\ 



+ — 

- 2-1- ^ 



H- 

— 2 — 

H—o J r 

02 a 3 

^a 2 

o 3 

di d^y 

a 3 ai 

^a 3 

Ol 

02 

aP J 


Baa — d j - 00000 1 
Bijij = ZCoiOiCojQj 


(4.14b) 

(4.14c) 


6 ofcOfc — 


2 (^didjdfz CLiQjjQjfo d^djd]^ didj’d ^ didjdfa H - ‘Zdi’dj(Ljc\ 

Lodfe \ ' 

x ^2idj / djd]£ d%djdfc dj,djdfa d%djdfc d'ldjdfc ~ l - 2 d'ldjdj^^ ^ 

/7? 

= 3, i ^ j ^ k. 


K D 

n 2 n 2- Dl J i j' 

u i a j 


(4.14d) 


Thus we have used three different definition of BR. The first one defined in ® 
imposes some restriction on the metric functions, namely for the BI spacetime is 
coincides with vacuum solution of Einstein equations. The second definition removes 
this restriction, but since BI metric admits non-trivial Ricci tensor, the BR in this 
case is not totally symmetric. Finally we gave the definition used by Bergqvist and 
Senovilla. It satisfies all the properties of BR and it is totally symmetric. In what 
follows we study the DEP and DSEP for BR in BI Universe. 
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5. Cosmological singularity and the dominant energy condition 

Recalling that a timelike geodesic is a world line for a particle moving without 
acceleration in the proper reference system we define the following: 

A spacetime is nonsingular if any timelike geodesics, or null geodesics, can be 
continued into the past and the future without bound, i.e., to infinite proper length 
for the timelike geodesics and to an infinite value of an affine parameter for the 
null geodesics. Such a spacetime is termed ’’causally, geodesically complete”. The 
requirements on the completeness are the minimum necessary so that the spacetime 
does not contain a singularity. It is necessary to point out that a spacetime not 
satisfying these requirements, however, one with a singularity, does not necessarily 
contain points with infinite curvature or with small hole. 

From physical point of view, of course, one ought to take as singular any space- 
time in which the geodesic world line of a particle cannot be continued without 
bound with respect to the proper time of this particle, for such a singular spacetime 
would lead to a violation of conservation laws. 

As applied to the cosmological problem, the Hawking-Penrose theorem reads as 
follows E3 

Theorem. A space-time M cannot satisfy causal geodesic completeness if the 
GTR (General Theory of Relativity) equations hold and if the following conditions 
are fulfilled: 

(i) The space-time M does not contain closed time-like curves. 

(ii) The energy condition (DEP) is satisfied at every point. 

The energy condition may be expressed as 

t a t a = 1 implies R a pt a t & < 0. (5-1) 

With Einstein’s equations 

Rap ~^9apR nTap 

El becomes 

t a t a = 1 implies T a pt a t^ > ^Tjf. (5.2) 

If, in an eigentetrad of T^, e denotes the energy density and pi, P 2 , P 3 denote the 
three principal pressure, then 15.21 can be written as 

£ + > 0 ; 

a 

£+Pa> 0, a = l,2,3. 

The weak energy condition is 

l a l a = 0 implies Rap^t 13 < 0, 


which is a consequence of ED- 
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(in) On each time-like or null geodesic 7 , there is at least one point for which 

K [a R b]cd[e K f] K c K d ± 0, (5.3) 

where K a is the tangent to the curve 7 at the given point and where the brackets 
on the subscripts imply antisymmetrization. If 7 is timelike, we can rewrite as 


R abcd K c K d ± 0. 

(iv) The space-time M contains either (a) a trapped surface, ( b ) a point P for 
which the convergence of all the null geodesics through P changes sign somewhere 
to the past of P, or (c) a compact space-like hypersurface. 

The DEP for the BI metric can be written in the form: 


Tg > Tla 2 + T 2 b 2 + T3 c 2 , 

(5.4a) 

T 0 ° > Tla 2 , 

(5.4b) 

T 0 ° > T| 6 2 , 

(5.4c) 

Tq > T|c 2 . 

(5.4d) 


In analogy with the Hawking-Penrose theorem stated above, Senovilla and oth¬ 
ers introduced DEP for the higher dimensional tensors. A detailed description of 
singularity theorems and their consequences can be found in the review paper by 
Senovilla “l 

The DSEP as defined by Senovilla reads: El 

Theorem. A rank-s tensor is said to satisfy the DSEP if 

T^X 1 -K‘ >0 (5.5) 

for any future-pointing causal vectors 

To justify its name the dominant DSEP obeys the following Lemma. 

Lemma. If a tensor T Ml ... Als satisfies the DSEP, then 


T 0...0 > llju-./J, V/zi,..., Us = 0,..., 71 — 1 (5.6) 

in any orthonormal basis {eT, |. 

It was also established in “ that any tensor satisfying the DEP possesses the 
following property: 

Property. satisfies DEP if and only if 


> 0 


(5.7) 


L ^1 .../Xs*' 

for any set If 1 ■ ■■l l f s of future-pointing null vectors. 

Let us now back to the Trank BR tensor and to check the DSEP for it. Since 
BR defined as 1)4.1 011 is a completely symmetric, trace-free 4-rank tensor, then it 
satisfy the DEP -. Therefore from the foregoing theorem, lemma and property for 
Bijki we can write: 


BabcdkikX^i > 0 , 
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B 0 ...o > \B abcd \, Va, b,c,d= 0, 1 , 2 ,3 (5.8) 

and 

BabcdFAlllt > 0 

Apparently, in view of 14.14bL Eq. 15.811 imposes some restriction on the metric 
functions, e.g., af < 1. But this is not the case, since we have used a coordinate 
basis to compute the BR tensor components I4.14all - I4.14dll and Eq. 15.01) refers to 
an orthonormal basis. We recall that in an orthonormal basis the components of 
the Weyl tensor obey the following relation: 

CoiOi = Cjkjk, i,j, k = 1,2,3 and i ^ j ^ k. (5-9) 

In view of (|5.9D the expressions (|4.14al) . (14.14b[) . (14.14cD and (|4.14dl) now read 

-^0000 = jkjk^ + ^ Ckiki ^ 4 “ ijij ^ (5.10a) 

Baa = -£?oooo> (5.10b) 

B%jij — ZCoiOiCojQj (5.10c) 

Bokok = Bijij. i,j, k = 1,2,3, ^ / j 7 ^ (5.10d) 

Thus in connection with the above Lemma relative to DSEP, Eq. 15.81) is fulfilled 

without any restrictions on the metric functions. 

6 . Conclusions 

In view of the importance of the BI model in the study of the present day Uni¬ 
verse we considered the most simple model with a perfect fluid as a source field. 
The corresponding solutions to the Einstein equations have been obtained. Three 
alternative definitions of Bel-Robinson tensor are considered. It is shown that the 
definition used by Deser et. al. is consistent with the Einstein equations when the 
source field is given by a vacuum only. The second definition used by Teyssandier 
is free from this restriction, but BR defined in this way is not totally symmetric. 
Definition used by Senovilla and Bergqvist does not suffer from this shortcomings, 
i.e., it has all the symmetries and the DSEP is satisfied. 
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